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ABSTRACT 

A sample of 18286 radio-loud AGN is presented, constructed by combining the 7th 
data release of the Sloan Digital Sky Survey with the NRAO VLA Sky Survey (NVSS) 
and the Faint Images of the Radio Sky at Twenty centimetres (FIRST) survey. Us- 
ing this sample, the differences between 'high-excitation' (or 'quasar-mode'; hereafter 
HERG) and 'low-excitation' ('radio-mode'; LERG) radio galaxies are investigated. A 
primary difference between the two radio source classes is the distinct nature of the 
Eddington-scaled accretion rate onto their central black holes: HERGs typically have 
accretion rates between one per cent and ten per cent of their Eddington rate, whereas 
LERGs predominatly accrete at a rate below one per cent Eddington. This is consis- 
tent with models whereby the population dichotomy is caused by a switch between 
radiatively efficient and radiatively inefficient accretion modes at low accretion rates. 
Local radio luminosity functions are derived separately for the two populations, for 
the first time, showing that although LERGs dominate at low radio luminosity and 
HERGs begin to take over at Li.4GHz ~ 10 26 WHz _1 , examples of both classes are 
found at all radio luminosities. Using the U/U max test it is shown that the two pop- 
ulations show differential cosmic evolution at fixed radio luminosity: HERGs evolve 
strongly at all radio luminosities, while LERGs show weak or no evolution. This sug- 
gests that the luminosity-dependence of the evolution previously seen in the radio 
luminosity function is driven, at least in part, by the changing relative contributions 
of these two populations with luminosity. The host galaxies of the radio sources are 
also distinct: HERGs are typically of lower stellar mass, with lower black hole masses, 
bluer colours, lower concentration indices, and less pronounced 4000A breaks indicat- 
ing younger stellar populations. Even if samples are matched in radio luminosity and 
stellar and black hole masses, significant differences still remain between the accretion 
rates, stellar populations, and structural properties of the host galaxies of the two 
radio source classes. These results offer strong support to the developing picture of 
radio-loud AGN in which HERGs are fuelled at high rates through radiatively-efficient 
standard accretion disks by cold gas, perhaps brought in through mergers and interac- 
tions, while LERGs are fuelled via radiatively inefficient flows at low accretion rates. 
In this picture, the gas supplying the LERGs is frequently associated with the hot X- 
ray haloes surrounding massive galaxies, groups and clusters, as part of a radio- AGN 
feedback loop. 

Key words: galaxies: active — radio continuum: galaxies — galaxies: jets — black 
hole physics — accretion, accretion discs 



1 INTRODUCTION 



Active Galactic Nuclei (AGN) are associated with the accre- 
tion of material onto supermassive black holes, of roughly a 
* Email: pnb@roe.ac.uk million to a billion solar masses, located near the centres of 
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essentially all massive galaxies (e.g. iMagorrian et al.lfl998l ). 
with a mass that correlate s strongly with the stellar mass 
(e.g. iMarconi fc Hundl2003l: iHaring fc Rixll2004r) or velocity 
dispe rsion (e.g. iGebhardt et al.l |2000| ; iFerrarese fc Merritd 
120001 ) of the surrounding galaxy bulge. It is now widely ac- 
cepted that the build-up of these supermassive black holes 
and that of their host spheroids are intimately linked. Ev- 
idence is growing that AGN activity may play an impor- 
tant role in the evolution of the host galaxy, with AGN out- 
flows being responsible for c ontrolling or terminat ing star- 
formation (e.g. see review bv lCattaneo et al"1l2009l ). 

AGN activity occurs in at least two different modes, 
each of which may have an associated, yet different, feed- 
back effect upon the host galaxy. The most commonly- 
considered mode of AGN activity is the 'standard' accre- 
tion mode associated with quasars. In this mode, which has 
been variously referred to as 'quasar-mode', 'cold-mode', 'ra- 
diative mode', 'fast-accretor', 'high-excitation', or 'strong- 
lined', material is accreted onto the black hole through a 
radiatively-efnci ent, optically-thick, geom etrically thin ac- 
cretion disk fe.g. IShakura fc Sunvaevlll973f ). These AGN ra- 
diate acro ss a very broad ra nge of the electromagnetic spec- 
trum (e.g. lElvis et al.|[l994 ) although a dusty structure sur- 
rounding the black hole and accretion disk, often referred to 
as a torus, obscures the em ission at some wa velengths when 
the AGN is seen edge-on (|Antonuccilll993l , and references 
therein). They are often associ ated with star-formation ac- 
tivity in the host galaxies (e.g. lKauffmann et al.|[2003al ). al- 
though possibly with del ays between the s t ar formation and 
the A GN activity (e.g. IWild et al.l l20ld ; iTadhunter et all 
l201ll . and references therein). A fraction of these AGN are 
radio-loud, possessing powerful radio jets that can extend 
for tens or hundreds of kpc. 

This radiatively-efncient accretion mode may be im- 
portant in curtailing star formation at high redshifts 
and setting up the tight relationship between black hole 
and bulge masses observed in the nearby Universe (e.g. 



and b ulge masses obser ved m tne nearby Universe (e.g. 
Silk & Rces 19981; lFabianlll999l : lKing||2003l ; iRobertson et al.l 



20061 ). Observational evidence that quasars can accelerate 



high-velocity winds is plentiful, although there is much 
debate as to whether these are thermal 'energy-driven' 
winds, or 'momentu m-driven' by radiation pressure (e.g. 
ICattaneo et al.ll2009l . and references therein). In radio- loud 
AGN, th e powerful radio jets may also shock-accelerate the 
gas (e.g. iBest et al.l|2000h . dri ving bipolar winds at speeds 
up to thousands of km/s (e.g. iNesvadba et al.|[20 08). Con- 
siderable uncertainty remains as to the mass and energy 
content of material that is driven out by these winds, the 
size-scale of the outflows in radio-quiet quasars, and the rel- 
ative importance of the different mechanisms that may drive 
the winds. 

There is a second mode of AGN activity, in which the 
accretion of material on to the black hole leads to little 
radiated energy, but can lead to the prod uction of highly- 
energe tic radio jets. It was first noted by iHine fc Longairl 
l|l979T ) that a population of low luminosity radio sources 
exist in which the strong emission lines normally found in 
powerful AGN were absent. It has since been shown that 
these radio sources also exhibit no accretion-related X-ray 
emission, nor infrared emission from a putative torus (e.g. 
lHardcastle et al1l2007l . and references therein), a nd are thus 
intrinsically different from the quasar-like AGN. IBest et all 



(2005af) showed that these low-luminosity radio sources are 
hosted by fundamentally different host galaxies to emission- 
line selected (quasar-like) AGN, in terms of their stellar 
masses and host galaxy properties. These AGN, which have 
been referred to as 'radiatively inefficient', 'radio-mode', 
'hot-mode', 'slow-accretor', 'low-excitation', or 'weak-lined', 
are believed to be fuelled through advection-dominated ac- 
cretion flows (ADAFs), which a re optically thin, ge ometri- 
cally thick, accretion flows (e.g. iNaravan fc YjHl995T ). They 
emit the bulk of their energy in kinet ic form through the 
radio jets (e.g. iMerloni fc Heinzl l2007i ); debate remains as 
to whether all of the energy of the jets is associated with 
accretion, or whether energy fr om the spin of the super mas- 
sive black hole is also tapped l|McNamara et al]|201ll . and 
references therein). 

Although the total cosmic contribution of the ener- 
getic output of jets is nearly two orders of magnitudes 
lower than that of rad i ation from the 'quasar-mode' AGN 
(e.g. ICattaneo fc Best! [2009), the jet energy is all de- 
posited locally to the system, potentially producing a very 
efficient feedback mechanism. This is most directly ob- 
served in the bubbles and cavities that radio-AGN are ob- 
served to evacuate in the hot hydrostatic gas haloes of 
their host galax i es or surrounding groups and clusters (e.g. 
Bohringer et al.l Il993l: ICarilli et al.lll994l ; iMcNamara et all 
2000l ; iFabian et al.ll2006l ). The energies estimated for the 
radio sources correlate well with the Bondi accretion rates 
expected from the hot gas (|Allen et al.|[200^ ); this suggests 
that this hot gas may form both the fuel for the radio source, 
and the repository of its energy, offering the potential for a 
feedback cycle. 

'Radio-mode' AGN have been widely used in galaxy 
formation models as a mechanism to switch off star forma- 
tion in the most massive galaxies, thus reproducing both the 
observed shape of the galaxy luminosity function and the 
"old, red and dead " nature of massiv e early-type galaxies 
i|Croton et al.ll2006l : iBower et al]|2006h . IBest et al l |2005bD 
showed that the prevalence of radio-AGN activity was a 
very strong function of the mass of the host galaxy, rising 
to over 30% in the most massive systems. Using a scaling 
relation between the r adio luminos i ty and the mechanical 
energ y of the jet (cf. Birzan et al. 2004; ICavagnol o et al.l 
l2010h . IBest et al.l l|2006h went on to show that the time- 
averaged energetic output of these sources is indeed suffi- 
cient to counter-balance gas cooling in early-type galaxies 
of all masses. On a larger-scale, radio-AGN are almost ubiq- 
uitous in the brightest cluste r galaxies of cool-core clusters 
l|Burnslll990l ; IBest et al.ll2007h , and have been invoked as the 
solution to both the 'cooling flow' and the 'entropy floor' 
problems in the intra-cluster medium of groups and clusters 
l|McNamara fc Nulsen|[2~007l . and references therein). 

The radiatively inefficient and radiatively efficient AGN 
clearly have fundamental differences, but the precise origin 
of these differences remains unclear. Some authors have ar- 
gued that it relates to the origin of the fuelling gas, with 
accretion of cold gas leading to a stable accretion disk 
and a radiatively efficient accretion, while the accretion of 
hot gas via the Bondi mechanism would produce the jet- 
domi nated radiatively inefficient AGN (e.g. lHardcastle et al.l 
l2007h . Others argue that the spin of the black hole is impor- 
tant {McNainara et al. 2 01ll ; lMartmez-Sansigre fc Rawlingd 
1201 ll ). A third hypothesis is that it is solely (or primarily) 
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driven by the Eddington-scaled accretion rate onto the black 
hole, with the ADAF mode occurring when the accretion 
rate is well below the Eddingt on limit. This was th e pre- 
diction in the original work of iNaravan fc Yil j 19951 ), and 
support for this picture has come from recent work indi- 
cating that broad-line AGN (ie. quasar-like AGN seen face- 
on) have lower limits t o their accretion rate s at around 1 
percent of E ddington (|Kollmeier et al. 2006; T rump et al.l 
2009a, 1201 lh . and indications that a switch between flat- 
spectrum radio quasars and BL Lac objects (which are be- 
lieved to be beamed versions of the radia tively- inefficient 
sources) also occurs a t that Eddington rate (|Wu et al.ll201ll ; 
iGhisellini et""aTll201ll l This hy pothesis is used by syn thesis 
models for AGN evolution (e.g. iMerloni fc Heindboosh that 
have been constructed based upon the two different accre- 
tion modes. 

A critical input to these AGN evolution models, and 
to understanding the evolvin g feedback role that AGN may 
play in galaxy evolution fcf. ICroton et a l. 2006; Bower et al.l 
[200g) is a full understanding of the different AGN popula- 
tions, their distribution in luminosity, their host galaxies, 
and their cosmic evolution. The cleanest method for select- 
ing samples of radiatively inefficient AGN is through radio 
selection, using the emission of their jets. Radiatively ineffi- 
cient AGN are detectable at other wavelengths: in particu- 
lar, around 30% of the population of 'X-ray bright optically 
normal' galaxies found in X-ray surveys show X-ray spectra 
with no absorptio n, yet no evidence of AGN activity at op- 
tical wavelengths (| Trump et a l. 2009b), and are interpreted 
as being radiatively inefficient AGN in which th e X-rays re- 
late t o a beamed component of the jet emission (Hart et al. 
2009). Nevertheless, only at radio wavelengths is the selec- 
tion function well-understood, and large samples can be con- 
structed. In addition, corresponding samples of radiatively 
efficient (radio-loud) AGN can be constructed at the same 
time in exactly the same manner, allowing direct compar- 
isons between the two. 

This paper presents a large sample of ra dio sources 
draw n from the Sloan Digital Sky Survey (SDSS; lYork et all 
2000), and compares the properties of radio-selected radia- 
tively efficient and inefficient AGN. The selection of the 
sample and the classification of the radio sources is de- 
scribed in SectionfJl Section investigates the nature of the 
radio sources: their accretion rates, their luminosity func- 
tion, and their cosmic evolution. Section [4] compares var- 
ious properties of the host galaxies of the two classes of 
sources. The results are discussed and conclusions drawn in 
Section [5] Throughout the paper, the cosmological parame- 
ters are assumed to have values of Q m = 0.3, J1a = 0.7, and 
H = TOkms^Mpc -1 . 



2 SAMPLE SELECTION AND PROPERTIES 

2.1 The overall radio source sample 

The sample of radio sources was constructed by com - 
bining the 7th data release (DR7; lAbazaiian et alj I2009T ) 
of the SDSS spectroscopic sample with the National 
Radio Astronomy Observatory (NRAO) Very Larg e Ar- 
ray (VLA) Sky Survey (NVSS: lOondon et al.1 \l9M ) and 
the Faint Images of the Radio Sky at Twenty centime- 



tres (FIRST) survey l|Becker et all 1 19951), broadly fol - 
lowing the techniques described by Best et all ((2005a) 
for the earlier SDSS DR2 sample. The parent sam- 
ple for the DR7 matching is the 927,552 galaxies in 
the value-added spectroscopic catalogues produced by the 
group from the Max Planck Institute for Astrophysics, 
and Johns Hopkins University (hereafter MPA-JHU), and 
available at http://ww w . mpa-garching . mpg.de/SDSS/ (cf. 
iBrinchmann et al. 2004). These galaxies were cross- matched 
with the NVSS and FIRST radio sources following the 
method of iBest et all ll2005al). b ut ad opting the improve- 
ment described bv lDonoso et ail ((2009) for identification of 
sources without FIRST counterparts. The cross-matching 
goes down to a flux density level of 5 mjy, which means that 
the sample probes down to radio luminosities of £i.4ghz ~ 
10 23 WHz~ 1 at redshift 2 = 0.1. The sample of detected 
radio sources is presented in Table [1] 

The next step was separation of the radio- AGN from 
star-forming galaxies. This has been improved since the DR2 
sample, and now makes use of an optimal combination of 
three different methods: the method based on 4000A break 
strengt hs and the r atio of radio luminosity to stellar mass, 
used in lBest et all l|2005ah : a method based on the ratio of 
ra dio to emission line lum inosity, similar to that presented 
in iKauffmann et all (l2008h: a standard 'BPT' emission-line 
diagno stic method (I Baldwin et all Il98ll : IKauffmann et ail 
l2003al ). Appendix IA1 provides full details of the combined 
method. The resultant classifications are provided in Ta- 
ble [1] The radio luminosity functions for star-forming galax- 
ies and radio-loud AGN separately provide broad confirma- 
tion of the success of the classifications (cf. Section 13. 1[) . 
Further tests have been carried out to ensure that none of 
the results of this paper is dependent upon the specific de- 
tails of the SF-AGN separation method. 

For the current paper, analysis is rest ricted to ra- 
dio s ources within the 'main galaxy sample' (|Strauss et all 
2002), comprising those galaxies with magnitudes in the 
range 14.5 < r < 17.77, and further restricted to the redshift 
range 0.01 < z < 0.3. Within this sample there are 9,168 ra- 
dio sources, of which 7,302 are classified as radio- AGN. The 
median redshift of the radio-AGN is z = 0.16, and 1,245 are 
located at redshift 2 ^ 0.1, to which redshift range some of 
the analyses are restricted. 

Properties of the radio source host galaxies are 
drawn from the value added catalogues of the MPA-JHU 
group. In particular, th ese include total stellar masses 
l|Kauffmann et alj|2003bh . accurate emission line fluxes, af- 
ter subtraction of the modelled stellar contin uum to account 
for u nderlying stellar absorption features (|Tremonti et all 
|2004| ). parameters determined directly from the spectra such 
as 4000A break strength s and galaxy velocity dispersions (cf. 
IBrinchmann et al. 1 120041 1 . and a compendium of basic pa- 
rameters from the imaging data such as galaxy magnitudes, 
colou rs, sizes and structural parameters (see lYork et all 
120001 . for more details). As noted by the MPA-JHU group, 
the formal line flux uncertainties quoted in the MPA-JHU 
DR7 catalogue significantly underestimate the true values 
(as determined by comparing derived line fluxes of sources 
observed multiple times), and so the line flux uncertainties 
have been scaled by the factors recommended by the MPA- 
JHU team. 
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Table 1. Properties of the 18286 SDSS radio galaxies. Only the first 20 sources are listed here: the full table is available electronically. 
The first three columns give the identification of the targetted galaxies through their SDSS plate and fibre IDs and the date of the 
observations. Columns 4 to 6 give the RA, Dec and redshift of the galaxies. Column 7 gives the integrated flux density of the source 
as measured using the NVSS. Column 8 provides the radio classification of the source, following Best et al (2005): class 1 are single- 
component NVSS sources with a single FIRST match; class 2 are single-component NVSS sources resolved into multiple components by 
FIRST; class 3 are single-component NVSS sources without a FIRST counterpart; class 4 sources are those which have multiple NVSS 
components. Where a galaxy has a central FIRST component, the integrated flux density and offset from the optical galaxy of that 
central FIRST component are given in columns 9 and 10. Column 11 provides a flag classifying the source as either a radio-loud AGN 
(1) or a star-forming galaxy (0), according to the criteria described in AppcndixfA] Column 12 indicates whether the source is included 
in the full statistical "main sample" studied in this paper (ie. SDSS main sample target with 14/5 ^ r ^ 17.77). Columns 13 and 14 
indicate whether selected radio-loud AGN are classified as LERGs or HERGs, respectively (if such classification is possible; sources with 
for both cases are unclassifiable using current data). 



Plate 


Julian 


Fibre 




RA 


Dec 


z 


i S TvT "\ T Q CI 
*-TN V DO 


Radio 


'-'FIRST 


Offset 


AGN 


Main 


LERG 


HEI 


ID 


Date 


ID 




(J2000) 




1.4 GHz 


Class 


1.4 GHz 






Samp 














(hr) 


(deg) 




(Jv) 




fly) 


(") 
\ J 










266 


51602 


5 


9 


784226 


-0.81043 


0.4486 


0.0069 


i 


0.0042 


1.76 


1 











266 


51602 


26 


9 


797071 


-0.34230 


0.1348 


0.0963 


i 


0.1010 


1.26 


1 


1 


1 





266 


51602 


100 


9 


742905 


-0.74164 


0.2038 


0.0068 


i 


0.0025 


0.46 


1 


1 


1 





266 


51602 


109 


9 


782474 


-0.25218 


0.1304 


0.0075 


i 


0.0043 


0.51 


1 


1 








266 


51602 


134 


9 


720425 


-0.54706 


0.3679 


0.0091 


2 


0.0000 




1 











266 


51602 


150 


9 


758252 


-0.36839 


0.0530 


0.0104 


1 


0.0010 


2.02 





1 








266 


51602 


179 


9 


746165 


-0.50828 


0.3693 


0.0059 


1 


0.0053 


0.16 


1 











266 


51602 


235 


9 


706746 


-0.00139 


0.1459 


0.0054 


1 


0.0049 


0.43 





1 








266 


51602 


439 


9 


720012 


0.41417 


0.0252 


0.0081 


3 


0.0000 
















266 


51602 


504 


9 


764199 


0.63871 


0.0303 


0.0052 


1 


0.0028 


2.09 





1 








266 


51602 


507 


9 


758241 


0.25554 


0.1291 


0.0275 


1 


0.0269 


0.32 


1 





1 





266 


51602 


550 


9 


776251 


0.46721 


0.4505 


0.0088 


1 


0.0058 


0.44 


1 











266 


51602 


554 


9 


787119 


0.66564 


0.0201 


0.0180 


1 


0.0131 


0.79 


1 


1 


1 





266 


51602 


559 


9 


786605 


0.70274 


0.0305 


0.0063 


1 


0.0045 


2.72 





1 








266 


51602 


577 


9 


785432 


0.73798 


0.2616 


0.0489 


1 


0.0094 


0.78 


1 


1 


1 





266 


51602 


617 


9 


805367 


0.78802 


0.2112 


0.0082 


1 


0.0078 


0.66 


1 


1 








266 


51630 


361 


9 


706815 


1.14969 


0.4498 


0.0523 


1 


0.0454 


0.56 


1 





1 





266 


51630 


529 


9 


772283 


1.08112 


0.5768 


0.0146 


1 


0.0112 


0.41 


1 











267 


51608 


19 


9 


907105 


-0.92869 


0.3583 


0.1848 


1 


0.1839 


0.77 


1 





1 





267 


51608 


34 


9 


944658 


-0.02334 


0.1391 


0.1660 


4 


0.0022 


0.27 


1 


1 


1 






2.2 High/low excitation classification of the 
sample 

A key requirement of the current analysis is an abil- 
ity to separate the radio source sample into the two 
fundamentally-different AGN classes. For consistency with 
previous works on radio source samples, the nomencla- 
ture of "high-excitation" and "low-excitation" radio sources 
(HERGs and LERGs) will be adopted in this paper to de- 
note these. 

The work of lLaing et al.l (| 19941 ). based on very pow- 
erful radio- AGN from the 3CR sample, suggested a fairly 
pronounced division between the two classes. Laing et al. 
classified as HERGs those sources which had the line flux 
ratio [OIII] 5007 / Ha > 0.2 and an equivalent width of 
the [OIII] line EWjoiii] > 3A (where em ission line equiv- 
alent widths are assigned positive values). iTadhunter et al] 
( 1998) similarly found that low-excitation sources in the 2Jy 
radio sample, with EWroni] < 10A, stand out from other 
galaxies in both their [OIII] 5007 / [Oil] 3727 line ratio 
and the ratio of emission line to radio luminosity. Lower lu- 
minosity radio samples, however, show a less-clear division 
between the two classes, since the emission line luminosi- 



ties of the high-e xcitation sources correlate s trongly with 
radio luminosity l|Rawlings fc Saunders! Il99ll ) and so be- 
come much weaker (and less easy to d istinguish from LERG 
lines) in low-luminosit y systems (cf. IZirbel fc Baumlll995l ; 
iKauffmann et alj|2008h . 

More recent analyses separating HERGs and LERGS 
have be en carried out for large samples of radio galax- 
ies (e.g. IButtiglione et al]|201fj ICid Fernandes et al]|2010l ; 
iBaldi fc Capettil2010h r largely based on emission line diag- 



nostics to separate Seyfert from LINER galaxies devised by 
iKewlev et al. 1 20061) . Note that a growing literature of work 
(e.g. ICid Fernandes et all 1201 ll . and references therein) in- 
dicates that a significant proportion of LINERs are not due 
to AGN activity, but rather the emission lines are photo- 
ionised by post-asymptotic giant branch (post-AGB) stars in 
old galaxies: such cases can be identifie d as having an equiv- 
alent width of the Ha line below 3A l|Cid Fernandes et al] 
l201ll ). In the radio galaxies studied here, the radio activity 
confirms that an AGN must be present, but this does not 
rule out the possibility that the weak emission lines could 
still have a post-AGB origin. 

IButtiglione et al] (|2010l ) have defined an 
"excitation index" parameter, combining four 
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Figure 1. The distribution of the HERG/LERG-classified radio 
sources on the [ QUI] equivalent width versus excitation index 
llButtiglione et al"1l2010T) plane, for galaxies with both parameters 
measured. This demonstrates the broad consistency of the two 
main approaches used to classify the radio sources. 



emission line ratios: EI = log 10 ([OIII]/H/3) — 
| [log 10 ([NII]/Ha) + log 10 ([SII]/Ha) +log 10 ([OI]/Ha)]. 
They demonstrate this parameter to be bimodal and use 
it to classify the galaxies, dividing the LERG and HERG 
populations at a value of EI — 0.95. For the SDSS radio 
galaxy sample defined in Section \2. 11 in many cases the full 
set of emission lines needed to classify the host galaxies 
via the Excitation Index are either not detected, or the 
signal-to-noise of the detections is low. Therefore a multiple 
approach was adopted to carry out the classifications, 
working down the following series of possibilities for each 
source until a classification was derived. 

(i) If all six emission lines were detected and the excita- 
tion index was at least la away from 0.95, the radio source 
was classified using the excitation index [783 LERG; 95 
HERG]. 

(ii) If four lines wer e relia bly detected for one of the in- 
dividual iKewlev et all l|2006() diagnostic diagrams and the 
source lay at least la from the division line, then that di- 
agnostic diagram was used for classification [330 LERG; 2 
HERG]. 

(iii) If the equivalent width of the [OIII] emission line was 
at least la above 5A then the source was classified as high 
excitation [65 HERG]. 

(iv) Classification options (i)-(iii) were repeated but with 
the la criterion removed [305 LERG; 30 HERG]. 

(v) If [Nil] and Ha emission line measurements were 
available for the source, t hen the [Nil] /Ha vs [OIII ] /Ha 
emission line diagnostic of ICid Fernandes et al.l (|2010T ) was 
used, using either a detection or a limit (if definitive) for the 
[OIII] emission line [769 LERG; 24 HERG] 

To illustrate the consistency of these cuts, Figure [1] 
shows the distribution of the classified galaxies on the 
EW[oiii] vs EI plane, for galaxies with both parameters 
measured. It is clear that the radio sources do show a good 
correlation and that both methods can provide suitable clas- 
sifications. There is a small population of LERGs with rel- 
atively high [OIII] equivalent widths but low excitation in- 



dices, and therefore a risk that classification method (iii) 
may lead to contamination of the HERG sample by such 
sources. However, such contamination is expected to be 
small: of the 65 sources classified by method (iii), 53 have 
excitation index measurements within la of 0.95, and all 
of the remainder have a limiting [OIII] /Ha > 1, strongly 
suggesting they are indeed HERGs. 

Despite using all of these different mechanisms, only 
about a third of the radio sources were able to be classi- 
fied. Many are simply undetected in emission lines, suggest- 
ing that they are likely to be LERGs - although in some 
cases this may just be due to the relative faintness of the 
source. To investigate a possible mechanism for classifying 
these sources, the left panel of Figure [2] shows the distribu- 
tion of [OIII] line luminosity versus radio luminosity for the 
classified sources. As discussed above, the two populations 
do occupy different regions of this planfl although there is 
significant overlap. The solid line in the figure represents an 
approximate lower limit of the distribution of the HERGs, 
and so the 3883 sources with emission line luminosities, or 
limits, below this can therefore be fairly securely classified 
as LERGs. As demonstrated in the right panel of Figure [2l 
this allows robust classification of all but four of the z < 0.1 
subsample of radio galaxies. 



3 FUNDAMENTAL PROPERTIES OF THE 
RADIO SOURCES 

3.1 Local radio luminosity functions of HERGs 
and LERGS 

Radio luminosity func tions were calculated in the s tandard 
way, as p = y~). 1/Vi |Schmidtlll968l ; rCondonlll989l ). where 
Vi is the volume within which source i could be detected. 
This is calculated as Vi = Vmax — V m i n , where Vmax and V m i n 
are the volumes enclosed within the observed sky area out 
to the upper and lower redshift limits, respectively, at which 
each source would be included in the sample. Redshift lim- 
its were determined by the joint radio and optical selection 
criteria, namely a radio cut-off of 5mJy and optical cut- 
offs of 14.5 < r < 17.77, as well as any imposed redshift 
limit for the analysis (e.g. 0.01 < z < 0.3). The sky area of 
the overlapping region between the SDSS DR7 spectroscopic 
survey and the FIRST radio survey, after removal of noisy 
regions around very powerful radio sources, was calculated 
to be 2.17 steradians. The summed radio luminosity func- 
tion of all radio sources, together with its separation into 
star-forming galaxies and radio-loud AGN, are provided in 
Table f2] and shown in Figure [3] an d are in excellent agree- 
ment with previous dete r minations jlMachalski fc Godlowsk: 
|2000l : ISadler et ai1l2002l : iBest et al.ll2005al : iMauch fc Sadie: 
2007). Uncertainties quoted are the statistical Poissonian er- 
rors only; at some luminosities these are so small that they 
will be under-estimates, with systematic errors dominating. 

Local radio luminosity functions were derived individ- 
ually for LERGs and HERGs; these are also tabulated in 
Table [2] and are shown in Figure [4] In order to minimise the 
potential influence of unclassified sources (typically found 



ki 

er 



1 Radio-quiet quasa rs would lie a bove and to the left of the 
HERG sources (e.g. lXu etahll 19991) . 



© 2011 RAS, MNRAS 000,[TJfl6] 



6 P. N. Best and T. M. Heckman 




Figure 2. Left: the [OIII] emission line luminosity versus radio luminosity for HERG/LERG-classified radio sources. The solid line 
indicates an approximate lower limit to the distribution of the HERGs, below which unclassified sources can be classified as LERGs with 
reasonable confidence. Note that only 20% (randomly selected) of LERGs are plotted to avoid over-crowding of the figure. Right: the 
same plot, but only for radio sources with z < 0.1, and also including the unclassified galaxies. It can be seen that all but four of the 
unclassified galaxies in this redshift range can be robustly classified as LERGs using the criterion derived from the left panel. 



Table 2. The local radio luminosity functions at 1.4 GHz, derived separately for the HERG and LERG populations. The first column 
shows the range of 1.4 GHz radio luminosities considered in each bin. The second and third columns show the total number of radio 
sources and the space density of these, in units of number per logio-L per Mpc 3 , detected out to z = 0.3. Columns 4 to 7 show the 
radio sources split into star-forming galaxies and radio-loud AGN. The eighth column gives the maximum redshift considered for the 
LERG/HERG analysis, in order to minimise the number of unclassified sources. The numbers and space densities of LERGs, HERGs, 
and unclassified sources, respectively, are given in columns 9 to 14. Uncertainties are statistical Poissonian uncertainties only. Note that 
the unclassified sources have negligible contribution compared to the LERG population, but for some bins of luminosity they would make 
a significant additional contribution if added to the HERG sample (see also Figure _J. 



logii.4GHz All radio sources Star-forming Radio-AGN LERGs HERGs Unclassified 

W Hz" 1 N log 10 p N log 10 p N log 10 P z max N log 10 p N log 10 p N log 10 p 



22.0-22.3 


297 


qo+0.03 


284 


o 1 1+0.03 
-o.ll_ .o 4 


13 


_ 4 44+0-15 
-0.23 


0.10 


12 


-4 46+ - 16 
4 -* D -0.25 


1 


_ 5 .74+o.2S 






22.3-22.6 


385 


-3 49+ 02 
-0.03 


357 


-3 <54+ 03 
°- d -0.03 


28 


_ 4 4,5+0.05 


0.10 


24 


_ 4 .1+0.05 


1 


-5 64+° 19 

°' D *-0.33 






22.6-22.9 


532 


o 27+0.02 
°-°'-0.02 


388 


-4 00+ ' 02 
^• uu -0.03 


144 


_ 4 45+O.04 
^•*°-0.04 


0.10 


138 


-4 46+ - 04 
4 - 4D -0.05 


5 


-6 01+° 16 


1 


-6.72 


22.9-23.2 


674 


-4 22+ ' 02 
-0.02 


298 


_ 4 .7+0.03 
^•°'-0.03 


376 


_ 4 40+0.02 
^•*°-0.03 


0.10 


339 


4 4Q+O.02 
^•* a -0.03 


9 


-6 09+ 13 






23.2-23.5 


882 


-4 56+ ' 02 

^• dD -0.02 


221 


r- -,0+0.03 
J - lo -0.03 


661 


-4 69+ ' 02 

-0.02 


0.10 


248 


, 7 q+0.03 

^•'°-0.03 


12 


-6 09+° 11 
°- UM -0.15 


3 


-6.66 


23.5-23.8 


1358 


_ 4 75+0.01 
^•'°-o.oi 


126 


-5 69+ 04 
°- oy -0.05 


1232 


-4 80+ ' 01 
^• ou -o.oi 


0.13 


377 


_ 4 qi+0.02 
-0.02 


12 


-6 42+ 11 
°'^ z -0.15 


10 


-6.48 


23.8-24.1 


1615 




52 


-6.35_™« 


1563 


4 01+O.02 
-0.02 


0.15 


454 


4 O i+0.07 
-0.09 


10 


-6 69+° 12 
°' DM -0.17 


8 


-6.77 


24.1-24.4 


1327 


-5.Q8±g 


19 




1308 


-5 09+ 01 

d - ua -0.02 


0.17 


427 


-5 19+ 002 

"-^ -0.02 


15 


-6 62+° 10 
°- DZ -0.14 


14 


-6.70 


24.4-24.7 


949 


r- 9 r- + 0.02 


3 




946 


5 26+ ' 02 
°- ZD -0.02 


0.17 


275 


r- oo + 0. 03 
-O.OO_ 03 


15 


-6 61+"' 10 
D ' D1 -0.13 


1 


-7.85 


24.7-25.0 


561 


5 K4+0.02 
°-°^-0.02 







561 


5 K4+0.02 


0.20 


228 


5 66+ 03 
-O.OD_ 03 


16 


-6 70+° 10 
,u -0.13 






25.0-25.3 


303 


5 ao+0.03 
°- OZ -0.03 







303 


5 ao+0-03 
°- OZ -0.03 


0.25 


206 


-5 91+ 03 

J - al -0.04 


21 


-6.76_°;°° 


5 


-7.63 


25.3-25.6 


103 


fi o 9 +0.05 
°-° -0.06 







103 


fi o 9 +0.05 
°- oz -0.06 


0.25 


57 


-6.441^7 


13 


-6.98_° ;°? 






25.6-25.9 


17 


fi r-o + 0.07 

-b.D»_ . og 







17 


fi r-o+0.07 
-O.o8_o.o8 


0.30 


29 


-6.82_°;°g 


17 


-6 95+° 10 
D ' yo -0.14 


1 


-8.55 


25.9-26.2 


12 


7 -1O + 0.12 







12 


-7.18_ Q 1? 


0.30 


9 


-7-4i±S:i 


3 


-7 64+ 17 

' -"^-0.28 






26.2-26.5 


3 


7 7S +0.21 
-'•'»_0.43 







3 


7 7S +0.21 
-'•' 8 -0.43 


0.30 







3 


7 7S +0.21 
-'■' 8 -0.43 







at higher redshifts) and yet retain sufficient volume for the 
rarer luminous sources, the upper limit of the redshift range 
used to calculate the radio source space density was in- 
creased with increasing radio luminosity, as indicated in Ta- 
ble__ Figure [4] illustrates the maximum effect that the resid- 
ual unclassified sources could have on the HERG luminosity 
function, even if all of them were HERGs: although there are 
small changes to some data points, the overall interpretation 
and conclusions are unaffected. The potential influence of 
unclassified sources on the radio luminosity function of the 
LERGs is completely negligible. 



Figure _] shows that, as expected, LERGs dominate the 
radio source population at relatively low radio luminosities, 
while the HERGs begin to dominate at the highest luminosi- 
ties, beyond Pi.4GHz ~ 10 26 WHz _1 . However, what is clear 
(and goes against standard simplifying assumptions in the 
literature) is that both populations are found across the full 
range of radio luminosities studied: even at radio luminosi- 
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Figure 3. The local radio luminosity function at 1.4 GHz de- 
rived separately for radio-loud AGN and star-forming galaxies. 
Filled points connected by solid lines indicate the data derived 
in this paper. For comparison, the results of Machalski & God- 
lowski (2000) using the Las Campanas Redshift Survey, Sadler 
et al. (2002) using the 2-degree field Galaxy Redshift Survey (2dF- 
GRS), and Mauch et al. (2007) using the 6-degree field Galaxy 
Survey (6dFGS) are shown. 
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Figure 4. The local radio luminosity function at 1.4 GHz, derived 
separately for the HERG and LERG populations. Note how both 
populations are found across the full range of radio luminosities 
studied. 



Taking account of the lower and upper redshift limits (see 
Section [3Tj) . the value of (V — Vmin) / (V max — Vm m ) was calcu- 
lated for each source. These were then averaged for HERGs 
and LERGs separately within radio luminosity bins; once 
again, to minimise the effect of unclassified sources, a vary- 
ing upper redshift limit was imposed for the analysis in each 
radio luminosity bin. The mean values are presented in Ta- 
bleland displayed in Figure [5] 

The HERG population displays clear evidence for cos- 
mic evolution at all radio luminosities studied, in the sense of 
there being a tendency for the sources to be located at larger 
distances than the mediar0. In contrast, the LERG popula- 
tion is broadly consistent with a mean V/Vm&x — 0.5 (with 
the possible exception of the highest radio luminosity bin). 
The LERGs therefore show little or no evidence for any cos- 
mic evolution. Combined with the results of Section fc.ll this 
result has important implications for understanding the evo- 
lution of the radio luminosity function as a whole, whereby 
the differential cosmic evolut ion seen betwe en powerful and 
less powerful radio sources (|Longairl 1 19661 ) may be driven 
by the switch in the dominant population with radio lumi- 
nosity (cf. Figure [3]). At high radio luminosities the HERGs 
dominate and, being a strongly-evolving population, lead 
to strong evolution of the overall radio source space den- 
sity at these luminosities (a factor ^t housand increase in 
space density out to redshift 2-3; cf. iDunlop fc Peacockl 
ll99d ; [Rigbv et aI1l201ll . and references therein). In contrast, 
at low radio luminosities the LERGs dominate the popula- 
tion, leading to the weak cosmic evolution seen in the low- 
luminosity radio populatio n (a factor 1.5-2 in c rease in space 
densi ty out toz~ 0.5, e.g. I Sadler et al. l l2007l ; lDonoso et all 
2009). Indeed, even this evolution may be driven to a large 
extent by the evolution of the HERGs: locally these con- 
tribute < 10% of the population at low luminosities, but if 
they increase in space density by an order of magnitude out 
to z « 0.5 (as they do at higher radio luminosities) then 
they would become comparable in numbers to the LERGs 
and lead to the observed doubling of the overall space den- 
sity. 



3.3 Eddington-scaled accretion rates of HERGs 
and LERGS 

As discussed in the introduction, a popular hypothesis for 
the difference between LERGs and HERGs relates to the 



ties around Pi.4GHz ~ lO^WHz" 1 the HERGs constitute a 
few percent of the overall radio-loud AGN population^- 

3.2 The cosmic evolution of HERGs and LERGS 

The cosmic evolution of the HERG and LERG populations 
can be individually investigated by using the V/Vmax test. 

2 Note that although there may be some contamination of 
the low-luminosity HERG population by radio-quiet quasars or 
Seyfert galaxies, the bulk of these are expected to be genuine 
radio-loud AGN: as detailed in Appendix A, the separation of 
star-forming galaxies from AGN was designed to exclude the 
radio-quiet quasars and Seyfert galaxies from the AGN category. 



3 Many HERGs arc classical-double radio sources, sometimes 
with no detected radio core. The possibility exists that sources 
without cores might be missed by the cross-matching procedure 
at low redshifts due to the large angular separation of the com- 
ponents. This would mean that 2 m j n should be higher than esti- 
mated, and could lead to an upward bias in (V — V m i n )/(V r max — 
Vmin)- In order to ensure that this is not the origin of the ob- 
served evolution of the HERGs, all HERGs with no detected ra- 
dio core were examined to determine the lowest redshift at which 
the cross-matching procedure would have included them in the 
sample, under the worst-case scenario that no additional FIRST 
components would be detected. For only three sources would this 
lead to a change in the estimate of z m ; n and these each lead 
to a change in the relevant (V — V m i n )/(V max — Vmin) value of 
significantly less than 0.01, which is negligible compared to the 
associated errors. 
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Table 3. V/V max determinations for the LERG and HERG pop- 
ulations separately, in various ranges of radio luminosity. In each 
luminosity range, an upper redshift limit (z max ) is defined for the 
analysis (and given in column 2) in order to reduce the number 
of unclassified sources to negligible levels. 
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Figure 5. The cosmic evolution, as a function of radio luminosity, 
of the HERG and LERG populations separately, as demonstrated 
using the V/V max test. 



Eddington-scaled accretion rates on to the black hole. The 
Eddington-scaled accretion rate can be estimated for the ra- 
dio sources by comparing the total energetic output of the 
black hole, calculated as the sum of the radiative luminos- 
ity and the jet mechanical luminosity, with the Eddington 
luminosity. 

The bolometric radiative luminosity of each radio 
source was estimated from the observed luminosity of the 
[QUI] 5007 emi s sion l ine, using the relation determined by 
iHeckman etafl (|2004T ): L rad = 3500L O iii- Where the [OIII] 
line was not detected, an upper limit to the radiative lu- 
minosity was set instead. The uncertainty on individual 
estimates of the bolometric radiative luminosity is ~ 0.4 
dex, from the scatte r around the L ra( j vs Loin relation 
(|Heckman et alj|2004h . 

The jet mechanical luminosity was estimated 
from the 1.4 GHz radio l umin osity, using the re- 
lation of ICavagnolo et all (|2010h . L mcch — 7.3 x 
10 36 (Li.4GHz/10 24 WHz- 1 )°' 70 W. This relation was 
determined using the energies associated with cavities 
evacuated by radio sources in the hot X-ray gas haloes 
of giant ellipticals, groups and clusters of galaxies. It is 
in broad agreement with minimum- energy synchrotron 
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Figure 6. Top: The distribution of Eddington-scaled accretion 
rates for the LERG and HERG populations separately. Analysis 
is limited to z < 0.1, in which redshift range essentially all radio 
sources could be robustly classified. For the LERGs, the solid line 
shows the best-estimate distribution, using the calculated values 
of the radiative luminosity. The dashed line shows the distribu- 
tion if, for all sources which are undetected in [OIII] , the 3<r upper 
limit to the [OIII] luminosity is used to calculate the radiative lu- 
minosity. The dot-dash line shows the distribution if the emission 
lines of all galaxies with Ha equivalent width below 3A are as- 
sumed not to have an AGN origin. These distributions therefore 
represent the allowed extremes of the LERG distribution, and 
show that the result is broadly similar in all three cases. For the 
HERGs, the dotted line is plotted with two different normalisa- 
tions: the upper line shows the fraction of sources relative to the 
total number of HERGs, while the lower line shows the fraction 
relative to the total number of LERGs to allow direct compari- 
son of numbers with the LERGs. Bottom: The distributions are 
shown considering only the radiative luminosity (ie. ignoring the 
mechanical energy in the radio jets). For LERGs, the dashed line 
indicates the maximal distribution calculated using either mea- 
surements of, or 3<t upper limit to, the [OIII] luminosity. The solid 
line indicates the subset of these which are measurements rather 
than limits. 



estimates of IWillott et alJ (|l999h . The scatter around the 
L m ech vs £i,4GHz relat ion is observed to be about 0.7 dex 
i|Cavagnolo et alj|20ld ). 

The black hole mass of the radio source host galaxy 
was estimated from the velocity dispersion of the galaxy 
(cr»), as measured in the SDSS spectrum, using the 
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well-established Mbh -o* relation. The determination of 
iTremaine et all |2002r i was adopted: log(AfBH/Afo) = 8.13+ 
4.021og(er*/200kms~ 1 ). The black hole masses thus derived 
define the Eddington limit for each radio source, L^dd = 
1.3 x 10 31 M B h/M q W. The intrinsic scatter in the Mbh-ct* 
relation is less than 0.3 dex. 

Combining these three relations, the Eddington-scaled 
accretion rate was derived for each radio source, as A = 
(irad + imcch) /^Edd- In cases where [OIII] was undetected 
(ie. signal-to-noise below 3) then the radiative luminosity 
was considered to be negligible compared to the mechani- 
cal luminosity, and the Eddington-scaled accretion rate was 
calculated from the mechanical luminosity alone. For these 
sources, a maximum value for the Eddington-scaled accre- 
tion rate was also calculated based on the upper limit to the 
[OIII] line luminosity (A max = (Z-rad,max + £ mec h) /-^Edd)- 
For sources where [OIII] was detected, but the equivalent 
width of the Ha emission line was below 3A, the possibil- 
ity was considered that the emission lines arise from photo- 
ionisation from post- AGB stars, instead of from the AGN 
(see Section 12.21 and ICid Fernandes et alj|201ll ). For these 
sources a minimum value for the Eddington-scaled accretion 
rate was determined, based only on the mechanical luminos- 
ity: Amin = Z/mcch /Z/Edd • 

The upper panel of Figure [6] shows the distribution of 
Eddington-scaled accretion rates for the HERG and LERG 
populations separately; analysis is limited to z < 0.1, in 
which redshift range essentially all radio sources are ro- 
bustly classified. For the LERGs, distributions are shown 
for each of A, A m i n and A max : the three distributions are 
similar, demonstrating that the high-A end of the distribu- 
tion of Eddington-scaled accretion rates for LERGs is ro- 
bustly described, regardless of the assumptions adopted (it 
should be noted that the low A end of the LERG distri- 
butions is influenced by the radio selection criteria which 
will lead to minimum detectable values of Z/ mcc h: the true 
distribution may extend to much lower values of A than indi- 
cated by Figure [6]). The Eddington-scaled accretion rates of 
the HERG and LERG populations are clearly fundamentally 
different: LERGs typically display accretion rates below one 
per cent of the Eddington rate, whereas HERGs typically 
accrete at higher Eddin gton rates (as is also found for radio- 
quiet AGN in SDSS; cf. IKauffmann fc Heckm"anll2009l ). This 
result very strongly implies that Eddington-scaled accretion 
rate on to the black hole is a primary factor in determining 
the nature of the accretion flow. The difference in the de- 
rived Eddington-scaled accretion rates between HERGs and 
LERGs is largely driven by the estimate of accretion rates 
onto the black holes (particularly the difference in Z, ra d, cal- 
culated from Loin) but also in part by HERGs typically 
being hosted by galaxies with less massive black holes (see 
Section 13]). 

For LERGs, the uncertainty in individual A measure- 
ments is ~ 0.7 dex, dominated by the intrinsic scatter in 
the Z/ m cch vs Li.4GHz relation. For HERGs, the radiative 
luminosity is generally much larger than the mechanical lu- 
minosity, and so the uncertainty in A measurements is ~ 0.4 
dex. These large uncertainties will make the observed dis- 
tributions of accretion rate broader than the intrinsic ones. 
Therefore, although there is some overlap in the accretion- 
rate distributions of the LERG and HERG populations, with 
examples of both classes being found between a few tenths 



and a few percent of Eddington, the results would be en- 
tirely consistent with a complete dichotomy whereby the 
LERG/HERG classification is enti rely determined by th e 
Eddington-scaled accretion rate (cf. iMerloni fc Hein3l2008l ) . 
Alternatively, the overlap in accretion rates between the two 
classes may imply that, although accretion rate is the pri- 
mary determinant, other factors may also be important in 
determining the nature of the accretion flow on to the black 
hole (black hole spin is an obvious candidate). The current 
dataset cannot definitively distinguish between these possi- 
bilities, because the lack of knowledge of the precise width 
and shape of the distribution of the uncertainties means that 
a robust deconvolution of the distributions is not possible. It 
is worthy of note, however, that if the observed LERG and 
HERG distributions are simply deconvolved with Gaussians 
of width 0.7 and 0.4 dex, respectively, then the two decon- 
volved distributions become entirely distinct, with a clean 
separation occurring at 1 percent of Eddington. 

Finally, the lower panel of Figure [5] shows the distribu- 
tions considering only the radiative luminosity of the sources 
(A ra d = Z/ ra( j/LEdd), i e - ignoring the mechanical luminosity 
of the radio jets. This is a somewhat unphysical comparison 
since it ignores the major energetic output of the LERGs, 
but is nevertheless included for comparison with other works 
in the literature which have defined the Eddington ratio in 
this way. The division between the two populations is even 
more clear in this analysis. 



4 THE NATURE OF THE RADIO SOURCE 
HOST GALAXIES 

Differences in host galaxy properties between the HERG 
and the LERG populations may provide interesting insight 
in to either the triggering mechanism of the radio activity, 
or the effect of the radio source on its host galaxy. Early 
work with relatively small samples compared radio galax- 
ies with strong emission lines against those with weak or 
absent lines, and suggested that the strong emission line 
sources were less luminous, with lower velocity dispersions, 
bluer colours an d lower mass-to-light ratios (ie. younger ste l- 
lar populations: [Smith fc Heck man 1989; S mith et al.| [l990). 
These results have be en confirmed wit h muc h larger ra- 
dio source samples by IKauffmann et ah I (|2008f l. who stud- 
ied the differences between SDSS-selected radio sources with 
and without detectable emission lines, and found the host 
galaxies of emission line radio sources to have lower stel- 
lar masses, lower velocity dispersions, lower 4000A break 
strengths, and stronger Balmer absorption features than 
those without emission lines. However, it is important to 
note that many LERGs do display emission lines, and so 
these results should be tested with properly defined LERG 
and HERG sa mples. 

iLin et all ([2010) also studied large samples of radio 
galaxies selected from SDSS, in their case separated into 
different radio morphological typ es (Fanaroff and Rile y 
Class 1 and 2 sources - FR1/2; iFanaroff fc Rilevl Il974h . 
These two radio morphological classes show a broad over- 
lap with the LERG and HERG classes, respectively, and 
therefore many differences between LERGs and HERGs 
may also be reflected by differences between FRls and 
FR2s. For this reason, many previous analyses of host 
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Figure 7. The distributions of stellar mass (upper left), black hole mass (upper right), 4000A break strength (lower left) and concentration 
index (C = -Rgo/R50> where Rgo and R50 are the radii containing 90% and 50% of the light, respectively; lower right), as a function of 
radio luminosity, for the LERG and HERG populations separately. The larger symbols indicate the mean values for each population in 
radio luminosity bins. 



galaxy properties, accretion rates, and luminosity func- 
tions have co ncentrated on dif f erence s between the FR1/2 
classes (e.g. | Ledlow fc Owenl Il99tj ; iGhisellini fe Celottil 
200ll ; ICao fc Rawlmdl2004lRigbv et al.ll200Sl ; lGendre et~ 
2010I . and references therein). However, there are substantial 
differences between the HERG/LERG and FR1/2 segrega- 
tions, since a sign ificant p opulation of LE RG FR2s exists 
|Laing et al.ll 1994) . Indeed. iLiii et all ||201Ch found the most 
significant differences when they compared the host galaxy 
properties of the most edge-brightened FR2s with strong 
emission lines against those of the other radio sources: this 
is closer to a HERG/LERG split. They found this FR2 sub- 
set to be host ed by lower mass galaxies, live in sparser envi- 
ronm ents (cf. iPrestage fc Peacocklll988l ; ISmith fc Heckmanl 
1990), and have higher accretion rates than the rest of the 
radio source population. Again this argues for the need to 
investigate clean LERG and HERG samples. 



Figure [7] shows the distributions of stellar mass, black 
hole mass, 4000A break strength, and concentration index 
(C — -R90/-R50, where -R90 and R50 are the radii contain- 
ing 90% and 50% of the light in the r-band) for the host 
galaxies of the LERG and HERG populations, as a func- 



tion of radio luminositjQ. It is immediately apparent that 
the HERG selection picks out host galaxies which are less 
massive and have lower black hole masses than those of the 
LE RGs, in lin e with the results of iKauffmann et all (|2008l ) 
and iLin et al.l (|2010T l. and earlier works. The 4000 A breaks 
strengths of the HERGs are also lower than those of the 
LERGs at all radio luminosities, indicating younger stellar 
populations (note that ther e is no issue of AGN li ght con- 
tamination; see discussion in lKauffmann et al.ll2008l ). At low 
radio luminosities the concentration indices of the HERGs 
are also lower. 

Many properties of galaxies in the local Universe corre- 
late strongly with the stellar mass. The difference in 4000A 
breaks strength between LERGs and HERGs in Figure [7] 
cannot therefore be properly interpreted until it is known 
whether it is simply driven by the lower typical stellar mass 
of the HERGs. To address this, matched samples of LERGs 
and HERGs were created. For each HERG, a search was 
made for LERGs which were matched to ±0.02 in z, ±0.1 
in logM, ±0.1 in logMan and ±0.25 in logl/Nvss- If at least 

4 Note that the same result for concentration index is found if 
an upper redshift limit of z = 0.15 is applied, indicating that any 
effects of seeing on the measurement of C are unimportant. 
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Figure 8. Histograms showing the distribution of radio source and host galaxy properties for samples of LERGs and HERGs which have 
been matched in redshift, stellar mass, black hole mass, and radio luminosity (with three LERGs matched to each HERG). The first 
four panels indicate the success of this matching, with the distributions of LERGs and HERGs agreeing well in these four parameters. 
The remaining panels show differences between the HERG and LERG populations in galaxy colours, sizes, concentration indices, 4000A 
break strengths, and Eddington-scaled accretion rates. Vertical dotted lines indicate the mean values for each population. 



three such LERGs were found, then three of these were ran- 
domly selected for the matched sample. If three matches 
could not be found then the HERG was excluded from the 
analysis. 

Figure [5] shows histograms of a variety of host galaxy 
properties for the matched LERG and HERG samples. The 
first four panels illustrate the distributions of redshift, stel- 
lar mass, black hole mass and radio luminosity for the 
matched samples, and demonstrate the success of the match- 
ing. The remaining five panels show the distributions of 
galaxy colour (g — r), size (R50 in kpc), concentration in- 
dex, 4000A break strength and Eddington-scaled accretion 
rate for the matched LERGs and HERGs. HERGs are seen 
to be bluer, smaller, less concentrated, and have lower 400fjA 
breaks than LERGs of the same stellar and black hole mass 
and radio luminosity. The offset in galaxy sizes is of marginal 
significance, as calculated by KS-tests on 1000 iterations 
of the random matching selection, but the offsets of the 
other four parameters are each significant at > 99.9% sig- 
nificance level. The difference in accretion rates is still ex- 
tremely strong in the matched samples (note that the me- 



dian of the LERG distribution is higher than in Figure [(J] 
because the matching with the HERGs pushes the LERG 
sample towards higher radio luminosities and lower black 
hole masses). The bluer colours and lower 4000 A breaks 
are consistent with the HERGs being associated with ongo- 
ing star formati on activity, as is also see n in similar radio- 
quiet AGN (e.g. IKauffmann et al.ll2003aT ). The smaller av- 
erage sizes and lower concentration indices may be related 
to the triggering mechanism of the sources. In particular, 
if LERGs are often fuelled by accretion from their hot gas 
haloes, then they might be more preferentially located at 
the centres of g roups or clusters (cf. IKauffmann et ai1l2008l ; 
iLin et al1 l2010). where host galaxies are typically larger and 
may have more extended (cD-type) light profiles. 



5 SUMMARY AND INTERPRETATION 

A large sample of radio sources drawn from the SDSS has 
been classified into high- and low-excitation radio galaxies, 
and the nature of these two different populations has been 
investigated. The main results are: 
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• Local radio luminosity functions have been derived sep- 
arately for the LERGs and HERGs for the first time. Both 
populations are found across the full range of radio lumi- 
nosities studied, although LERGs dominate the population 
at low radio luminosities and HERGs at high luminosities. 
The two populations appear to be switching in dominance 
at L1.4GHZ ~ lO^WHz" 1 . 

• HERGs show evidence for strong cosmic evolution at 
all radio luminosities, while LERGs are consistent with little 
or no evolution. This difference, coupled with the changing 
population mix as a function of radio luminosity, helps to 
drive the strong luminosity-dependence of the evolution seen 
in the overall radio luminosity function. 

• The accretion rates of the HERGs and LERGs are fun- 
damentally different. HERGs typically have accretion rates 
between one and ten percent of Eddington. LERGs, in con- 
trast, predominantly accrete at rates of below one percent 
of Eddington. 

• HERGs are hosted by galaxies of lower mass and lower 
black hole mass than LERGs of the same radio luminos- 
ity. The host galaxies of HERGs are also bluer and have 
lower 4000A break strengths than LERGs of the same mass 
and radio luminosity, indicating the presence of associated 
star formation. LERGs are larger and have more extended 
light profiles than HERGs, consistent with the interpretation 
that they are more likely to be hosted by central galaxies of 
groups and clusters. 

These results are consistent with the developing picture 
of radio-loud AGN, in which HERGs are fuelled at relatively 
high rates in radiatively-efficient standard accretion disks by 
cold gas, perhaps brought in through mergers and interac- 
tions, and with some of the cold gas leading to associated 
star formation. The requirement for significant cold gas sup- 
plies means that these sources are much more prevalent at 
earlier cosmic epochs where merger rates and gas fractions 
were larger. In contrast, LERGs are fuelled at relative low 
rates, through radiatively-inefficient accretion flows, largely 
by gas associated with the hot X-ray haloes surrounding 
the galaxy or its group or cluster (although gas from any 
other source fuelling the black hole at low accretion rates 
would also lead to a LERG). Regardless of whether the ac- 
cretion occurs directly from the hot gas through the Bondi 
mechanism, or after the hot gas has cooled, this gas source 
allows the setting up of an AGN feedback cycle, since this 
is the gas direct ly affected by any radio-AGN activity (cf. 
iBest et al1l2006l ). In this picture, LERGs will be associated 
with massive galaxies, which have old passive stellar pop- 
ulations. Given that massive galaxies show little evolution 
out to z ~ 1, the lack of cosmic evolution seen in the LERG 
population is as expected. 

One of the most interesting results coming out of this 
study is the almost distinct nature of the accretion rate 
properties of the LERG and HERG classes. This is in line 
with earlier indications co mparing FR1 and FR2 sources 
l|Ghisellini fc Celottil lioOlh , and with more recent results 
for a switch in accretion rate between B L Lacs and flat- 
spect rum radio quasars l|Wu et all l201ll ; iGhisellini et al.l 
2011). This result fits in well with theoretical calculations 
of accretion modes onto black holes (jNaravan fe Yilll995h . 
and is in line with what is found in galactic X-ray binary 
systems. X-ray binaries display three well-defined spectral 



states (e.g. iFender et al.l [20o3 . and references therein): a 
'low/hard' state in which the source is characterised by hard 
X-ray emission, low-power radio jets are ubi quitous, and 
the r adio and X-ray luminosities are correlated l|Gallo et alj 
I2003L possible analogue to the radiatively inefficient AGN 
state); a 'high/soft' state dominated by a thermal X-ray 
component characteristic of a standard thin accretion disk 
(possible analogue to radio-quiet quasar- like AGN); a short- 
lived 'intermediate' or 'transition' state between the two, 
where both the thermal accretion disk and powerful radio 
jets are seen (possible analogue to radio-loud quasar-like 
AGN). The switch between these different accretion states 
in X-ray binary systems has been shown to depend upon 
the accretion rate onto the black hole, with the switch oc- 
curri ng at between 1 and a few percent of the Eddington 
rate (M accaron el l2003l) - very similar to the value derived in 
this paper for radio-AGN . 

iKording et al.l l|2006l) demonstrated that spectral states 
analogous to those of X-ray binaries could be identified in 
local AGN, while other analogies between the black holes 
in X-ray binaries and those in AGN have been uncovered 
through examinations of the relationships between radio and 
X-ray luminosities and black hole mass (the so-called 'Fun- 
damental plane of black hole activity'; iMerloni et al.l [20031 ; 
iFalcke et al.1l2004 ). On the basis on those results, and the 
assumption that not only do AGN show these same three 
spectral states but that there is also the same dependence 
on accretion rate between them, synthesis models of AGN 
have been constructed (e.g. lMerloni fc H einz 200^). The re- 
sults in this paper offer strong evidence in support of these 
AGN synthesis models, and the broad principles of the anal- 
ogy of AGN with X-ray binary systems, by both confirm- 
ing previous indications of a lower limit to the accretion 
rate of radiatively efficient AGN, and demonstrating that 
an upper limit of around a percent Eddington applies to the 
radiatively-inefficient population. 
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APPENDIX A: SEPARATION OF 
STAR-FORMING GALAXIES AND 
RADIO-LOUD AGN 

This appendix provides details of the techniques used to 
classify the selected SDSS radio sources as either star- 
forming galaxies or radio-loud AGN. It is important to clar- 
ify that the attempt made here is to classify the origin of the 
radio emission, and hence that radio-quiet AGN are clas- 
sified together with the star-forming galaxies, rather than 
contaminating the radio-loud AGN category. 

Three mechanisms for separating the two classes of 
sources are considered: (i) using the relationship between 
the 4000A break str ength and the rat io of radio luminosity 
per stellar mass fcf. lBest et aLll2005aT ). hereafter referred to 
as the 'D4000 vs -^rad/M' method; (ii) using emission line 
diagnostics, in particular the ratio of [Qui] 5007 and H/3 
line f luxes, and that of [Nil] 6584 and Ha (cf. iBaldwin et al.l 
Il98ll ), hereafter referred to as the 'BPT' method; (iii) using 
the relation between the Ha emission line luminosity and 
the radio luminosity - the 'Lho vs L ra d method. The pre- 
cise divisions adopted for each of these three methods are 
described in Section I All and are illustrated in Figure IA1I 
Following this, the manner in which the results of these 
methods are combined to produce an overall classification 
is explained in Section \A2\ 
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Figure Al. The location of radio sources on the three classifica- 
tion lines used to separate radio-loud AGN from galaxies where 
the radio emission is powered by star-format ion. The top plot is 
the 'D4000 vs L ra d/yV method, developed bv lBest et al] ll2005al V 
The middle plot shows the widely-used 'BPT' emission line ra- 
tio diagnostic. The lower plot shows the relationship between Ha 
and radio luminosity. In all plots, the dotted lines indicate the 
division used for that classification method. Sources plotted as 
red diamonds are classified as radio-loud AGN in the overall clas- 
sification, while star-forming galaxies appear as black crosses. In 
the lower plot, arrows indicate upper limits to the Ha luminosity. 
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Table Al. The table presents the numbers of sources in each combination of classifications by the three different classification methods, 
together with the overall classification adopted and the rationale for that classification. It should be emphasised that 'AGN' in the 
classification refers to sources classified as radio-loud AGN, whereas 'SF' refers to sources where the radio emission is powered by 
star-formation (albeit that a radio-quiet AGN may also be present). 



D4000 vs 


BPT 




No of 


Overall 


Rationale for overall classification 


^rad/M 




vs L ra( j 


Sources 


class 




AGN 


AGN 


AGN 


1078 


AGN 


Unambiguous AGN 


AGN 


AGN 


?? 





- 


No sources 


AGN 


AGN 


SF 


847 


AGN 


Mostly close to cut line in Ljy a but clear AGN in other plots 


AGN 


?? 


AGN 


8053 


AGN 


All weak-lined, reliable AGN 


AGN 


?? 


?? 


1938 


AGN 


Weak-lined AGN. Majority high-z. hn a limit nearly allows AGN classification 


AGN 


?? 


SF 


222 


AGN 


Clear AGN in D4 oo- Close to Lj^a cu t; offset from bulk of SF population 


AGN 


SF 


AGN 


25 


AGN 


Radio excess. Near BPT cut line. Probably radio-loud AGN with associated SF 


AGN 


SF 


?? 





- 


No sources 


AGN 


SF 


SF 


40 


SF 


Close to D4000 cut line. Most have high extinction. 


?? 


AGN 


AGN 


42 


AGN 


Clear AGN; almost all have high L ra( j and high D4000 but no mass estimate. 


?? 


AGN 


?? 





- 


No sources 


?? 


AGN 


SF 


55 


SF 


Most likely scenario is radio-quiet AGN with associated SF 


?? 


?? 


AGN 


1166 


AGN 


Mostly z > 0.3. Almost all have high L ra( j and high -D4000 but no mass estimate. 


?? 


?? 


?? 


1767 


AGN 


All but six have z > 0.3 and L rad > 10 25 W^Hz _1 


?? 


?? 


SF 


12 


SF* 


Properties mostly appear like SF galaxies; see footnote 


?? 


SF 


AGN 







TxTr^ crtiirffic 


?? 


SF 


?? 







No sources 


?? 


SF 


SF 


29 


SF* 


Consistent classification, but see footnote 


SF 


AGN 


AGN 


37 


AGN 


Mostly high z, high L ra( j. Probably radio-loud AGN with SF activity 


SF 


AGN 


?? 







No sources 


SF 


AGN 


SF 


1582 


SF* 


Mostly low luminosity; almost certainly radio-quiet AGN 


SF 


?? 


AGN 


73 


AGN 


Mostly high z, high L ra( j. Probably radio-loud AGN with SF activity 


SF 


?? 


?? 


12 


SF* 


Clear SF in D4000 but too distant for emission line classifications 


SF 


?? 


SF 


35 


SF 


Consistent classification 


SF 


SF 


AGN 


19 


SF* 


Mostly low L rat i SF galaxies near Lh„ cut line, but see footnote. 


SF 


SF 


?? 







No sources 


SF 


SF 


SF 


1224 


SF 


Unambiguous SF 



* In these classes, there exist small populations of objects with L ra( j > 10 24,5 WHz 1 and z > 0.3. The properties of these subsets of 
high L ra d, high-z galaxies appear much more like radio-loud AGN. Therefore, although the bulk of these objects are classified overall as 
SF galaxies, the subsets of -L ra d > lO^WHz" 1 and z > 0.3 sources arc classified as radio-loud AGN. 



Al The three classification methods 

The 'P 4000 vs L ra d/M-' method was developed bv lBest et al.l 
(2005a|), on the basis that star- forming galaxies with a wide 
range of star formation histories occupy a similar locus in 
this plane. This is because, for star-forming galaxies, both 
L r ad/M and D4000 depend broadly on the specific star for- 
mation rate of the galaxy. Radio-loud AGN will have en- 
hanced values of L ra d and are thus separable on this plane. 
iBest et alJ (|2005aT ) adopted a division line of 0.225 higher 
in D4000 than the track produced by a galaxy with an ex- 
ponentially declining star formation rate of 3 Gyr e-folding 
time. They demonstrated that this method was generally 
very successful by comparison w ith other method s used 
in the literature. Further work by Kauffm ann et all (|2008l ) 
indicated that this selection line may be somewhat too 
shallow at low values of L ra( j/M, misclassifying some star- 
forming gal axies as AGN. Tak ing account of this, the divi- 
sion line of IBest et afl (|2005ah was modified such that for 
irad/M < 12.2, a straight-line cut is adopted with equation 
£>4ooo = 1-45 — 0.55(L rad /M — 12.2). Most radio sources 
are classifiable by this method, but w 17% of sources are 
not, due to the lack of a mass estimate in the MPA-JHU 
SDSS value-added catalogues. These unclassifiable sources 
are almost always the highest redshift sources, z > 0.3, se- 
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lected in samples other than the 'main galaxy sample' (ie. 
with r > 17.77), although small numbers of lower redshift 
sources where the automated algorithm had failed are also 
present. 

The 'BPT' method has been widely used by a number of 
authors to divide star-forming galaxies from AGN, since the 
different spectrum of the ionising r adiation leads to di fferent 
emission line ratios. As outlined by Best et al. (2005a), how- 
ever, there are problems relating to its direct use for identifi- 
cation of a clean radio-loud AGN sample. In particular, it is 
well-established that AGN and star-formatio n activity are 
closely related (e.g. Kau ffmann et al.l l2003al l. Radio-quiet 
AGN may therefore be detected in the radio on the basis 
of their on-going star-formation activity and identified as 
AGN on the basis of their emission line ratios, thus con- 
taminating the radio-loud sample. Nevertheless, the emis- 
sion line ratios can offer useful information in many cas es. 
Here, the division proposed bv lKauffmann et al.l (|2003ai ) is 
adopted, namely that galaxies with log([Oin]/H/3) > 1.3 + 
0.61 / (log([Nll]/Hoj) - 0.05) are classified as AGN. Just un- 
der 30% of the radio source sample are classifiable by this 
method, but the rest lack detections (or definitive limits) for 
at least one line. 

The 'Lho vs Lrad'' method is based on the premise that, 
for star-forming galaxies, both the Hcv luminosity and the 
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radio luminosity provide a direct measure of the star forma- 
tion rate and these properties are therefore correlated. For 
radio-loud AGN, a far higher radio luminosity to Ha lumi- 
nosity ratio is observed, and so the location of galaxies in the 
Lif e vs L ra d plane can be used to identify radio- loud AGN 
(cf . iKauffmann et ai1l2008h . Ideally the extinction-corrected 
Ha luminosity would be used, since dust attenuation can 
decrease the Ha-to-radio luminosity ratio of star-forming 
galaxies. However, this is only possible for the subset of 
galaxies for which H,9 is also detected, to allow a dust atten- 
uation estimate from the Balmer decrement, and even then 
the attenuation correction that should be used is different for 
star-forming galaxies and AGN. Here the observed Ha lu- 
minosity is therefore used instead, with an adopted division 
line of log{L Ha /L e ) = 1.12 x (logfL.^/WHz- 1 ) - 17.5). 
This cut is designed to be conservative, in the sense that 
galaxies classified as radio-loud AGN by this cut should be 
secure AGN, whereas the SF class may be contaminated by 
some AGN close to the cut line. Almost 80% of the radio 
sources are classifiable by this method, either using the ob- 
served Ha luminosity, or because the upper limit on that 
luminosity allowed clear classification. 

A2 Combination into a final classification 

For each of the three classification methods, every radio 
source is classified as a radio-loud AGN, classified as hav- 
ing its radio emission associated with star formation, or is 
unclassified. There are therefore 27 possible combinations 
of classification. Table [ATI indicates the number of galaxies 
in each of these 27 different classes and the resulting clas- 
sification adopted. These overall classifications are arrived 
at by examining in detail the properties of the galaxies in 
each class, especially where different classification methods 
disagree. For example, in most cases where objects miss the 
D4000 vs L ra d/M classification due to the lack of a stellar 
mass estimate, the values of D4000 and I/ ra( j are both so 
high that the galaxy will be classified as an AGN for any 
plausable value of the stellar mass and so that classification 
can be securely adopted. Similarly, in some classes with dis- 
puted classifications, if the sources lie very clearly within the 
AGN regime in two of the plots but just on the SF side of 
the cut line in the third, then that adds weight to an overall 
AGN classification. The rationale for the final classification 
is included in Table IAT1 

Figure [3] shows the radio luminosity functions de- 
rived for star-forming galaxies and radio-loud AGN. The 
smooth nature of the derived luminosity functions, es- 
pecially at low and high luminosities where small num- 
bers of misclassifications would have pronounced effects, 
offers broad support to the success of the classification 
scheme, as d oes the comparison with prev i ous determi- 



( 2002) and lMauch fc Sadler! |2007l ). The slight offset of the 
AGN line from previous determinations below a radio lumi- 
nosity of L fa 10 23 WHz _1 may be caused by the different 
classification of radio-quiet AGN, which here are classified 
along with the star-forming galaxies. In any case, objects of 
these low luminosities are excluded from most of the analy- 
ses in this paper. 



nations 
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